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Abstract. We discuss electron scattering in a one-dimensional delta barrier potential with 
either time-dependent coupling constant (classical model) or a coupling constant that is linear 
in a boson coordinate (quantum model). We find an exact continued fraction solution and Fano 
like resonances in the transmission coefficient. In the quantum model, energies for perfect 
transmission exist below the first sideband threshold. 



1. Introduction and Model 

Simple models for the interaction between fermions and bosons continue to be fascinating, 
as often very non-trivial results can be obtained from even the most primitive Hamiltonians. 
Exactly solvable models for the interaction of photons or phonons with electrons in quantum 
dots [[]]] or quasi-one-dimensional systems [[|,|3[] provide the best examples, as they often 
provide a deep insight into rich and complex physics. 

In this contribution, we re-consider a simple model for a single electron of mass m in 
one dimension that interacts with a delta-barrier through a coupling parameter that itself is a 
dynamical quantity. The Hamiltonian is written as 

H = ^— + 5(x){g + g^a) + a}) + Vta)a. (1) 
2m 1 " ' 

Here, a) creates a boson of frequency Vt and g\ [a) + a] is a dynamical contribution on top 
of the static coupling constant g . The constant zero point energy is omitted since it merely 
shifts the energy scale by Vt/2. The lattice version of this model was originally introduced 
by Gelfand, Schmitt-Rink and Levi [Q] years ago in the study of tunneling in presence of 
phonons, and was shortly afterwards compared to a corresponding time-dependent classical 
Hamiltonian [|J], the continuous version of which reads 

2 

H c i{t) = ^- + 6{x) {g + 20i cos(fit)} . (2) 
2m 

H c \ (t) is obtained as the interaction picture Hamiltonian of Eq.(JI[) with respect to Hb = Via) a, 
after replacing the boson operators by — a — 1. 

In its time-dependent version, Eq.(Q) has subsequently been used as a model for 
scattering in quasi-one-dimensional quantum wires by Bagwell and co-workers [[3[|5[], who 
found Fano-type resonances in the transmission coefficient as a function of the energy of 
an incident electron. It soon turned out that the scattering properties of this Hamiltonian are 
quite intriguing as they very much depend on the relative sign and strength of the two coupling 
parameters g and g\. The interplay between evanescent modes, quasibound states and the 
behaviour of the transmission amplitude in the complex energy plane [|§,[7p have been studied 
recently. 

Our focus here is on the quantum version Eq. (|T]) of the model and its peculiarities in 
comparison with H d . It turns out that beside transmission zeroes, there are points of perfect 
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transparency in the Fano resonance that only appear in the model H but not in H d . Perfect 
transmission and Fano resonances have been found recently in the transmission of phonons 
through non-linear chains without delta impurities |]||, §|. Although not discussed in detail 
here, these results indicate that there still is rich and non-trivial behavior to be discovered 
from models like Eq.(jl]). 



2. Transmission Coefficient 

The total wave function |^) of the coupled electron-boson system can be expanded in the 
oscillator basis {|n)} as 

oo 

(z|tt) = 5>nO*OM (3) 
n=0 

with wave function coefficients ip n (x) depending on the position x of the electron. We solve 
the stationary Schrodinger equation at total energy E > 0, implying a scattering condition for 
the electron part of the wave function in demanding that there is no electron incident from the 
right. For i^O, the ip n (x) are superpositions of plane waves if E is above the threshold for 
the n — th boson energy, 

ip n (x < 0) = a n e iKx + b n e~ iknX 

ip n (x > 0) = t n e iknX , k n = y/E - ntt, E > ntt, (4) 
whereas normalizabale evanescent modes occur if E is below the threshold, 

ip n (x < 0) = b n e KnX 

ip n (x > 0) = t n e~ KnX , n n ee y/nQ - E, E < nQ. (5) 

Here and in the following we set h = 2m = 1. We impose the condition that the boson is in 
its ground state for an electron incoming from the left, 

a n = <5n,o, (6) 

where we set the corresponding amplitude A = A to unity. Continuity of ^ n (x) at x = 
yields a n + b n = t n for all n, whereas the jump in derivative of ip n (x) across the delta barrier 
leads to a recursion relation for the transmission amplitudes t n , 



giVntn-i + (go - 2ij n )t n + gWn + lt n+1 = - 2i^ n 5 nfl (7) 
where the 7 n are real (imaginary) above (below) the boson energy nfl, 

In — k n 6(E — nf2) + iK n 9(nVt — E). (8) 
The total transmission coefficient T(E) is obtained from the sum over all propagating modes, 

[E/n] , , F] 

t(e)= E TmUm 2 , (9) 

n =o fc ol^J 

where the sum runs up to the largest n such that k n remains real. 



3. Matrix Representation and Continued Fractions 

Although Eq.@ is a finite sum, its evaluation requires the solution of the infinite recursion 
relation Eq.(Q) due to the fact that the propagating modes are coupled to all evanescent modes. 



M 



(10) 



Transmission Through a Quantum Dynamical Delta Barrier 

The transmission amplitudes can be determined from the linear equation 

Mt = a, t = (to,ti,t 2 ,...)> a = (-2i7o,0,0,...) 
f g - 2i7 \f\g x \ 

\f\g x go - 2z7i \f2gx 

\/2gi g - 2i>y 2 ' ' • 

V ••• J 

Numercally, this can easily been solved by truncation of the matrix M. Alternatively, one 
can solve Eq.(|i"0l) recursively which actually is numerically more efficient. For example, the 
result for the zero-channel transmission amplitude t (E) can be written in a very intuitive 
form: defining the 'Greens function' Gq{E) by 

G (E) = [-2ir /0 (E)+g ]- 1 , (11) 

one can write to(E) with the help of a recursively defined 'self energy' (E), 

t (F \- -2z7o(£) yWfrr\- m\ 

01 } Go\E)-W{Ey { } Go\E-Nn)-EW+V(Ey 

In fact, using ^ n {E) = ^(E — nCl), the self energy EW (E) can be represented as a continued 
fraction 

£«(£) = ^ — « . (13) 



G^(s-n) ^ 

G 1 (E-2fi) 



3^ 



4fl 2 



This demonstrates that t (E) depends on g\ only through g\. 

Truncating the matrix M to a N x iV matrix corresponds to the approximation that sets 
E^(£) = and recursively solves Eq. Q for E (Ar_1) (£?) down to EW(£). For example, 
truncating at = 2 we obtain the approximation 

An important observation can be made with respect to the stability of our theory for large 
coupling constants g\. In fact, the truncation at N + 1 is only consistent if the truncated self 
energy YP^'{E) is a small correction to the inverse 'free propagator', 

Ng\j\Gl\E - Ml) | < \G \E — (N — l)n)\, (15) 

which by use of Eq. ( |TT]) a t large N implies Ngf < ANil or g\ < 2y/Q. The tridiagonal 
form of the matrix, Eg. (|10|), actually implies that the method based on the recursion Eq. (0) 
is perturbative in the coupling g\ to the boson. We conjecture that for g\ above the critical 
value, the perturbation based on the oscillator basis {\n)} used here breaks down. A similar 
breakdown of numerical approaches that start from a weak coupling regime in single boson 
Hamiltonians is known from the Rabi Hamiltonian []T0|], i.e. the coupling of a single boson 
mode to a spin 1/2. 

We mention that the lattice version of the present model would be a natural starting point 
for a more detailed analysis of the strong coupling (small polaron) limit. 
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Figure 1. Transmission coefficient through a dynamical one-dimensional delta barrier with 
repulsive (go > 0, left) and attractive (go < 0, right) static part, cf. Eq. ([!]) and (||). E is the 
energy of the incident particle. 



4. Comparison to the Classical Case 



The recursion relation corresponding to Eq. ([7]) for the classical time-dependent Hamiltonian, 
Eq. (Q), was derived and discussed by Bagwell and Lake [Q], 

gitn-l + (90 ~ ^ln)tn + ftWl = ~ ^iJnfi, U = 0, ±1 , ±2, . . . (16) 



Here, t n is the coefficient of the time-dependent electron wave function in photon side-band 
n, where n runs through positive and negative integers n. In further contrast to the recursion 
relation Eq. (0), there are no factors ^fn and y/n + 1 multiplying the coupling constant g\. 
This latter fact is an important difference to the quantum case where these terms lead to the 
factors iV that multiply g\ in the self energies Y,( N \E), Eq. (|T2"|). This difference is eventually 
responsible for the breakdown of the perturbative approach for large g 1 in the quantum case. 

A continued fraction representation of t (E) for the classical case has been derived 
recently by Martinez and Reichl [Q. The corresponding matrix defining the transmission 



amplitudes t c 



t_ 2 , t-l, to, t\, t 



in the classical case is the infinite tridiagonal 



matrix M c \ with g — ij n on the diagonal and g\ on the lower and upper diagonals, 

( ■. n \ 





'•• #0-227-1 gi 
g x g - 2z7o 







gi 







gi 

go - 2«7i 



(17) 



V u •• ••/ 

Fig. (|l]) shows a comparison between the transmission coefficient T(E), Eq. @, for the 
quantum and the classical barrier. In the repulsive case with < g\ < go, the dynamical part 
of the barrier is essentially a weak perturbation to the unperturbed (#i = 0) case. Additional 
structures (cusps) appear at the boson (photo side-band) energies nQ although the overall 
T(E) -curve resembles the (g\ = 0) case. 

The more interesting case occurs for barriers with an attractive static part, g < (Fig. 
([]]), right). A Fano type resonance appears below the first threshold E = Vt where the 
transmission coefficient has a zero in both the classical and the quantum case. In the classical 
case, this is a well-known phenomenon [§]: the transmission zero for weak coupling (small 
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Figure 2. Logarithmic plot of reflection coefficient R = 1 — T for dynamical delta barrier 
with static repulsive (go > 0, left) and attractive (go < 0, right) core. 



gi) shows up when the Fano resonance condition 

2k 1 (£) + <7o = (18) 

is fulfilled.There, the energy of the electron in the first side channel (n = 1) coincides with 
the bound state of the attractive delta barrier potential, E — Vt = —g^/A. In the quantum case, 
the self energy in Eq.(|T2"|) diverges at the zeros of T(E), 



0. 



(19) 



For g x -> 0, YP^{E) -> T$ =2 (E) = gl/(2K X {E) + g ), cf. Eq.(0), and the two conditions 
Eq.Q) and Eq.([l9[) coincide. 



5. Perfect Transparency 

The most interesting feature in the scattering properties of the dynamical quantum barrier 
is the appearance of an energy close to the first channel (n = 1) threshold where perfect 
transmission T(E) = 1 occurs. This is clearly visible in the vanishing of the reflection 
coefficient, R(E) = 1 — T(E), in the logarithmic plot Fig. @. For a repulsive static part, 
go = 0.3, this occurs at an energy below the energy where the reflection coefficient comes 
close to unity, and above that energy if the static part is attractive (g = —0.9). In contrast, in 
the classical case the reflection coefficient never reaches zero in neither the repulsive nor the 
attractive case. This contrast becomes even more obvious in a two-dimensional plot where 
the zeros in R correspond to 'ridges' in the g -E plane, cf. Fig. @. 

Perfect transparency (R = 1 — T = 0) can be understood by considering the transmission 
amplitude t (E) which determines the total transmission below the first sideband threshold. 
Recalling that t (E) = —2ik / (—2ik + g — (£/)), in the quantum case the transmission 
coefficient becomes unity when 

g -^ 1 \E)=0. (20) 

Our exact continued fraction expression for the self energy, Eq.(JT3|), implies that for < E < 
fi, Tj^(E) is real because G () 1 (E — nil) = — E + g is real for n > 1. The condition 

Eq.(]20|) then means that the self energy exactly renormalizes the static part g of the scattering 
potential to zero. 

For small g x , we can use our perturbative expression corresponding to truncating the 
matrix M, Eq.([10|), to a two-by-two matrix. The perfect transparency condition Eq.(|2"0|) then 
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Figure 3. Density plot of In R (reflection coefficient) for the quantum delta barrier at g\ = 0.5. 
Exact solution from Eq. ( |l2| ) (left), from the N = 2 truncation Eq. (14) (center), and from the 
classical model Eq.(^j) (right). The light 'ridges' correspond to curves of perfect transmission 
T, cf. Eqs. (HI) and Eqs. (^ 



becomes 

9o- £\ , = 0, 0<E<n, (N = 2 truncation.), (21) 

2k 1 (E) +g 

which determines the position of the perfect transmission energy. The solution of Eq.(pTl) 
defines two curves in the E-g -p\ane with perfect transmission for < E < fi, 

g = -Vtt-E ± y/n-E + gl. (22) 

These two curves can be clearly identified in the logarithmic density plots of the reflection 
coefficient R = 1 — T, cf. Fig. @. The N = 2 approximation to the transmission amplitude, 
Eq.([B]), turns out to reproduce these features quite well even at moderate coupling constants 

9i- 

6. Conclusions 

The above analysis of the two models Eq. ([]]) and Eq. (|2|) has revealed some interesting 
differences between scattering properties of simple electron-boson models and their classical 
counter-part. The strong coupling limit of the quantum model and its extension to more 
complicated situations like multi-channel scattering remain to be explored. 

This work was supported by the UK project EPSRC GR44690/01, and the UK Quantum 
Circuits Network. 
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